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rich isotopes of oxygen and fluorine, and employ "bare" and secondary renormalized nucleon-nucleon 
interactions. For the nucleon-nucleon interaction from chiral effective field theory at order next- 
to-next-to-next-to leading order, we find that the coupled-cluster approximation including triples 
corrections binds nuclei within 0.4 MeV per nucleon compared to data. We employ interactions 
from a resolution-scale dependent similarity renormalization group transformations and assess the 
validity of power counting estimates in medium-mass nuclei. We flnd that the missing contributions 
due to three-nucleon forces are consistent with these estimates. For the unitary correlator model 
potential, we find a slow convergence with respect to increasing the size of the model space. For 
the G-matrix approach, we flnd a weak dependence of ground-state energies on the starting energy 
combined with a rather slow convergence with respect to increasing model spaces. We also analyze 
the center-of-mass problem and present a practical and efficient solution. 
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I. INTRODUCTION 

In the last decade, ab-initio nuclear structure calcu- 
lations have made great progress [i|-@|- Light nuclei up 
to carbon or so can now be described in terms of their 
nucleonic degrees of freedom, realistic nucleon-nucleon 
(NN) forces (i.e., those that fit the available body of NN 
data with a x^ « 1 per datum) augmented by a three- 
nucleon force (3NF). One of the major advances is due 
to the systematic construction of nuclear forces within 
chiral effective field theory (EFT) [To|-[l5|. Within this 
approach, the unknown short-range part of the nuclear 
force is systematically encoded in terms of low-energy 
constants and contact terms, while the long-range part 
of the interaction stems from pion exchange. One of the 
hallmarks of this approach is the "power counting", i.e., 
an expansion of the nuclear Lagrangian in terms of the 
momentum ratio Q/A^. Here, Q denotes the typical mo- 
mentum scale at which the nucleus is probed, while A,,, 
denotes the high-momentum cutoff scale that limits the 
applicability of the effective field theory. Within this ap- 
proach, three-nucleon forces appear naturally at next-to- 
next-to-leading (N^LO) order {Q/K^Y , and four-nucleon 
forces get introduced at N^LO, that is at order (Q/A^)** 
in terms of the momentum scale and the cutoff A^. 

The chiral interactions have been tested in very light 
systems of mass A=3, 4 through precise few-body calcu- 
lations [3, 01, and in p-shell nuclei within no-core shell 
model (NCSM) calculations [g, 0|- In these calculations, 
the NN interaction was taken up to N'^LO, while the 
three-nucleon interaction was taken at its leading order 
{Q/Axf (N^LO). Within the NCSM calculations, the 
"bare" chiral interaction needs to be renormalized, due 
to the size of the model space. Lattice calculations pro- 
vide a different approach for the implementation and so- 
lution of chiral interactions in the nuclear many-body 



problem Q. Presently, such calculations employ two- 
and three-nucleon forces at order N^LO, and they target 
p-shell nuclei. The coupled-cluster method is an alterna- 
tive approach that is particularly suited to study the sat- 
uration properties of chiral nucleon-nucleon interactions 
in medium-mass nuclei such as ^°'^^Ca and ^®Ni Q, o r 
even heavier nuclei. Using a Gamow basis expansion [l6| , 
this method can also include continuum effects and de- 
scribe rare isotopes that are in the vicinity of nuclei with 
closed (sub)shells [13, El- 

The "bare" chiral interaction models are still rela- 
tively "hard" interactions, that is the interaction has 
nonzero matrix elements even for high momentum trans- 
fers Q ^ kp, where kp ~ 1-35 fm~^ denotes the Fermi 
momentum at nuclear saturation. Using for example ma- 
trix elements of the chiral interaction by Machleidt and 
Entem [13| , we extend typically the momentum integra- 
tions up to 8 fm~^. For a model space consisting of os- 
cillator wave functions, a simple estimate shows that the 
number of oscillator shells required for the solution of a 
nucleus with radius R and an interaction with momen- 
tum cutoff A is about N « XR (see Sect. |TV] below for 
details). Recall that the number of single-particle states 
(single-particle j-shells) grows as N^ in the m-scheme 
(and as N"^ in a j-coupled spherical scheme). It is there- 
fore clear that high-momentum interactions require large 
model spaces. This makes wave-function-based solutions 
of the nuclear many-body problem a challenging task. 
Within the coupled-cluster approach, one can overcome 
this difficulty through a spherically symmetric reformu- 
lation of the method. This approach enables us to con- 
sider much increased model spaces, including even up to 
20 oscillator shells. The obvious alternative consists of 
lowering the momentum cutoff. 

The momentum cutoff of the interactions can be de- 
creased by a similarity renormalization group (SRG) 



transformation of the NN interaction [l9[, or by inte- 
grating out high-momentum modes within a renormal- 
izati qn g roup transformation or a similarity transforma- 
tion j2Ct l2ll | . Such transformations preserve all proper- 
ties (such as phase shifts and bound states) of the NN 
interactions up to the cutoff. The SRG transformations 
yield interactions that are band diagonal in momentum 
space, and the band width is characterized by a momen- 
tum scale. Eigenstates of these interactions are localized 
in momentum space, and the low-momentum states are 
thus decoupled from the high-momentum physics. The 
low-momentum interactions Viow k are characterized by a 
momentum cutoff beyond which all interaction matrix el- 
ements quickly approach zero. We denote the momentum 
scale of the SRG interactions and the cutoff of the T^ow fc 
interactions by A. Starting with A = A^, one can thus 
generate a one-parameter family of interactions from the 
chiral NN interaction models by performing SRG trans- 
formations, or by integrating out high-momentum modes. 
We will use this approach in order to study the scale- or 
cutoff dependence of the nuclear binding energy, and to 
examine aspects of the power counting. Throughout this 
work we make the underlying assumption that a com- 
plete description within chiral effective field theory would 
yield an accurate description of atomic nuclei. Thus, we 
attribute any scale dependence directly to the neglected 
three-nucleon forces and other neglected high-order terms 
of the interaction. 

The more traditional way of dealing with "hard" in- 
teractions consists of the computation of the so-called 
G-matrix !:22:i|, which is based on a Green's function 
approach [23] using normally unperturbed propagators. 
The G-matrix depends thus on the starting energy that 
is employed in its construction. Analytical arguments de- 
mand that this dependence weakens as increasingly larger 
model spaces are considered, and it vanishes for infinite 
spaces. Within the spherical coupled-cluster approach, 
we are able to investigate the convergence properties of 
the G-matrix that is constructed from chiral interaction 
models. 

Another approach to renormalize "hard" interactions 
is the unitary correlator method (UCOM). Here a unitary 
transformation is constructed to remove the hard core 
and short-range contributions of the tensor force by an 
appropriately formulated correlation operator [2J|. This 
method has seen several applications, but its convergence 
and saturation properties have not yet been studied in a 
framework that includes substantial wave function corre- 
lations. We fill this gap in this work. 

It is the purpose of this paper to employ "bare" 
and renormalized nucleon-nucleon interactions (employ- 
ing SRG and low-momentum techniques, the G-matrix 
and the UCOM method), to compare and analyze their 
convergence and saturation properties and their impact 
on nuclear structure. In addition to this task, this paper 
also contains significant supplemental information that 
could not be pre sented in several recent short communi- 
cations (sl. [l8l l25l - [27| . In particular, we present a detailed 



study and practical solution of the center-of-mass prob- 
lem, and study the evolution of single-particle energies in 
neutron-rich oxygen and fluorine isotopes. 

This paper is organized as follows. In Section|lTl we in- 
troduce spherical coupled-cluster theory. We address and 
resolve questions regarding the center-of-mass problem in 
Section Hm In Section ITVl we compute the binding ener- 
gies of the nuclei ^He, ^^O, '"^Ca, and ""^Ca starting from 
"bare" chiral NN interactions. Section |V] is dedicated 
to the evolution of single particle energies in neutron- 
rich oxygen isotopes. In Section IVTl we study ''"Ca for 
several momentum scales A of the SRG interaction, and 
examine the power counting. Section IVIII focuses on the 
starting-energy dependence and the convergence proper- 
ties of the G-matrix for ''He and ^^O. We analyze the 
convergence properties of the UCOM interaction in Sec- 
tion IVHll We finish with our Summary. Some technical 
details of the spherical coupled-cluster method are rele- 
gated to the Appendix. 



II. SPHERICAL COUPLED-CLUSTER THEORY 

In this section we g ive an outline of the coupled- 
cluster method |28l - [39| and introduce an angular mo- 
mentum coupled formulation of the coupled-cluster equa- 
tions. First we outline coupled-cluster theory for the 
computation of ground state energy of closed-shell nuclei 
within the so-called CCSD and A-CCSD(T) approaches. 
Thereafter, we introduce the equation-of-motion theory 
for the calculation of ground and excited states in closed- 
and open-shell nuclei. We introduce also the spherical 
formulation of coupled-cluster theory. Finally, we dis- 
cuss how to calculate expectation values of observables 
in coupled-cluster theory. 



A. Coupled-cluster theory for closed-shell nuclei 

Coupled-cluster theory is based on the similarity trans- 
formation 



H ^ e^^He^ 



(1) 



of the normal-ordered Hamiltonian H . Here, the Hamil- 
tonian is normal-ordered with respect to a product state 
|(/)o) which serves as a reference. Likewise, the particle- 
hole cluster operator 



T = Ti+T2 



Ta 



(2) 



is defined with respect to the reference state. The k- 
particle /c-hole (fcp-fch) cluster operator is 



Ti - ^ \^ ,ai...a. t „t „ fj. 



(fc! 



ii,...,ifc;ai,...,afc 



(3) 



Here and in the following, the indices i,j, k, . . . label oc- 
cupied single-particle orbitals while a,b,c, . . . label unoc- 
cupied orbitals. The most commonly used approxima- 
tion is coupled-cluster with singles-and-doubles excita- 
tions (CCSD) where T « Ti -h T2. The unknown am- 
plitudes tf and tfj in Eq. ^ are determined from the 
solution of the coupled-cluster equations 



= 
= 






(4) 
(5) 



Here \(t>f) 
state, and 



i,ab\ 



'0) is a Ip-lh excitation of the reference 
„ ) is a similarly defined 2p-2h excited state. 
The CCSD equations ^ thus demand that the reference 
state \(f>o) has no Ip-lh and no 2p-2h excitations, i.e., 
it is an eigenstate of the similarity-transformed Hamilto- 
nian ([ij in the space of all Ip-lh and 2p-2h excited states. 
Once the CCSD equations are solved, the ground-state 
energy is computed as 



E = 



\H\M 



(6) 



Coupled-cluster theory has the following virtues. First, 
the similarity-transformed Hamiltonian can be evaluated 
exactly since the Hausdorff-Baker-Campbell expansion 



H 



H + [H, T] 

1 
+ 3! 
{He-) 



[[i/,T],r] 



[[[ff,r],T],r] 



(7) 



terminates exactly at fourfold-nested commutators for 
two-body Hamiltonians. Here, the last term in parenthe- 
ses, (...)c, indicates that only those terms enter where 
a cluster operator T is connected to the Hamiltonian 
H to the right. Since every cluster operator is linked 
to the Hamiltonian, coupled-cluster theory by construc- 
tion fulfills Goldstone's linked cluster theorem and yields 
size-extensive results. This is particularly important in 
applications to medium-mass nuclei. Second, within the 
CCSD approximation, the computational effort scales as 
n^ni^j where Uo and n„ denote the number of occupied and 
unoccupied orbitals in the reference state \4>o), respec- 
tively. Thus, the computational effort is much smaller 
than within the configuration interaction approach (or 
shell-model in nuclear physics) which exhibits a factorial 
scaling as function of the chosen single-particle space and 
number of nucleons. 

Coupled-cluster theory with inclusion of full triples 
(CCSDT) [io'l is usually considered to be too compu- 
tationally expensive in most many-body systems of con- 
siderable size. Therefore triples corrections are usually 
taken into account perturbatively using the non-iterative 
CCSD(T) approach described in Ref. ^41,]. Recently, 
a more sophisticated way of including the full triples 
known as the A-CCSp(T) approach, has been developed 
by Taube et al. ^^. In the A-CCSD(T) approach 
the left-eigenvector solution of the CCSD similarity- 
transformed Hamiltonian is utilized in the calculation of 



a non-iterative triples correction to the coupled-cluster 
ground state energy. The left eigenvalue problem is given 
by 

{i;\KH ^ E{i^\K , (8) 

were A denotes the de-excitation cluster operator 

A = 1 + Al + A2 , (9) 

with 

Al = ^AXat, (10) 



A2 = - ^ X'lt,abaaa\a] . 



(11) 



ij^a^b 



The unknowns, A^ and A^-'j,, result from the ground-state 



solution of the left eigenvalue problem 



They are 



utilized together with the cluster amplitudes, tf and tf^, 
to compute the energy correction due to triples clusters 
as 



AK, = 



1 



(3 



n2 



J2 {m{Fhp + v)N\i^tjk] 



ijkabc 



X 7^A^t^k\iyNT,)c\ 



(12) 



-ijk 



Here, F^p denotes the part of the normal-ordered one- 
body Hamiltonian that annihilates particles and creates 
holes, while 



^abc — 

'ijk 



— In > Jjj t Jkk Jaa Jbb J a 



(13) 



is expressed in terms of the diagonal matrix elements 
of the normal-ordered one-body Hamiltonian F. In the 
case of Hartree-Fock orbitals, the one-body part of the 
Hamiltonian is diagonal and F^p vanishes. The subscript 
C denotes the connected part of the operator, and \ipfj'^) 
is a 3p-3h excitation of the reference state. 



B. Equation-of- motion coupled-cluster theory for 

ground and excited states of closed- and open-shell 

nuclei 

The CCSD and the A-CCSD(T) approaches are known 
as single-reference coupled-cluster methods (SR-CCM), 
and therefore they are particularly well suited for nu- 
clei with an expected closed (sub-) shell structure. By 
adding or removing particles to a closed-shell nucleus, 
we move into regions of nuclei with open-shell struc- 
ture, and clearly the SR-CCM methods are not well 
suited for the description of these nuclei. Most nuclei 
are of open-shell character and in order to study and 
predict properties like the evolution of shell-structure 
as we move towards the drip line in various isotopic 
chains, we clearly need to go beyond the SR-CCM class 
of approaches. There exists a variety of coupled-cluster 



methods that have been specifically designed to address 
the structure of open-shell systems, see Ref. 33] for 
an overview of some of these methods. Most methods 
are either based on a mult i- reference formulation of the 
coupled-cluster method (MR-CCM) (see Ref. ^M] for a 
discussion of various SR-CCM and MR-CCM) ^_or an ex- 
tension of equation-of-motion (EOM) theory [4^ based 
on the coupled-cluster method (EOM-CCM) (see, e.g., 
Refs. Hi, 113). The beauty of the EOM-CCM is that 
it has the simplicity and transparency of SR-CCM such 
as the CCSD and the A-CCSD(T) approaches discussed 
above, and that the method allows for systematic im- 
provements and extensions. 

In order to extend our ab initio coupled-cluster pro- 
gram beyond closed (sub-) shell nuclei, we use the EOM- 
CCM approach. We give a brief outline of the method 
in the following. The basic idea behind EOM [4a| is to 
calculate states \tpk) of the nucleus B by acting with 
an excitation operator J7fc on the ground state \ipo) of 
a closed-shell reference nucleus A. If ftk is a particle 
number-conserving operator (B — A), Qk will generate 
excited states of the nucleus A. In the case of fifc not 
conserving the number of particles {B ^ A), fl^ can gen- 
erate ground and excited states of the open-shell nucleus 
B with respect to the reference nucleus A. In EOM-CC 
theory, the ground state of the closed-shell reference nu- 
cleus A is given by the coupled-cluster wave function, 
\ipo) — exp{T)\ipo). In case we wish to calculate excited 
states of the closed-shell nucleus A {B = A), we define 
the excitation operator Qk 



■Rfc' 



Rt = ro + }^ rtala, + t ^ 



Y^ 



4 ^1^ r^^alalaju^- 



(14) 



and by truncating at the two-particle-two-hole excitation 
level, we get the standard EOM-CCSD approach. In case 
we wish to calculate ground and excited states of a nu- 
cleus with a nucleon added or removed from a closed-shell 
nucleus {B = AzLl), we define the excitation operators 
^k = Rt^^ with, 

<"'^ = E ^^^^ + J E ^H^^^i + • ■ • ' (15) 



and truncating at the two-particle-one-hole and two- 
hole-one-particle levels, we get the the particle-attached 
and particle-removed equation-of-motion coupled-cluster 
methods (PA/PR-EOM-CCM), see, for example, Refs. 
m 113. The PA-EOM-CCM adds a particle to a nu- 
cleus with A nucleons by creating Ip, 2p-lh, 3p-2h, . . ., 
excitations on the ground state of the nucleus with A 
nucleons. Similarly, PR-EOM-CCM removes a particle 
from A by creating Ih, lp-2h, 2p-3h, . . ., excitations on 
the ground state of A. By multiplying the equation for 
the nucleus A, H\(j)o) = Eo\(j)o), from the left with ilk, 



and subtracting it from the equation for the nucleus B, 
Hflk\(t>o) = Ek^k\4>o)i we obtain the following equation. 



[H,nk] |0o) - (i/r!fe)p|0o) =^fcf^fc|0o) 



(16) 



with H = 



,-T 



He^ being the similarity-transformed 



Hamiltonian and uik = Ek — Eq the energy of the nu- 
cleus B relative to the ground state of the nucleus A. 
Eq. (|16p defines a right eigenvalue problem for the exci- 
tation amplitudes in Eqs. p4)) and ([15]), and is usually 
solved using iterative eigenvalue algorithms. 



C. Spherical coupled-cluster theory 

For "hard" interactions such as the "bare" N'^LO in- 
teraction, wave function methods need to employ very 
large model spaces in order to yield converged results. 
For spherical reference states (nuclei with closed major 
shells or closed subshells) , it is therefore better to employ 
the spherical symmetry to further reduce the number of 
unknowns, that is, the number of cluster amplitudes. For 
such nuclei, the cluster operator of Eq. (I2|) is a scalar un- 
der rotation, and depends only on reduced amplitudes. 
Thus, 



Ti 



E*i«x«..)^°' 



(17) 



and 



T, 



hJjiajbJ 



*ia;('^)(<x<)^'^-(«.. x«..)^'^- (18) 



Here, we employed the usual notation for spherical ten- 
sors, and ji and ja denote the spin of the occupied and 
unoccupied subshells, respectively. It is clear that the 
similarity-transformed Hamiltonian is also a scalar un- 
der rotation, and it is straightforward to work out the 
CCSD equations within this formulation. Details are 
given in the Appendix. The computational cost of PA- 
EOM-CCSD is like that of CCSD, i.e., ntnl, where n„ is 
the number of unoccupied orbitals and Uo is the number 
of occupied orbitals. In order to reduce the memory and 
computational cost related to the basis size, we write the 
excitation amplitudes Rj. as spherical tensors of rank 
J and projectionM, 



ja 

+1 E '^r('^^"^)fH«4 



(19) 



jjjajbJab 



Uj 



{J) 

M 



Ri^-'\j,M)^Y.''"^J)^nm.5'j^ 
ji 

+1 E ^iu,(^."^)K®^^ 

jijjjaJij 



(20) 



iJ) 

M 



Here, we wish to solve for the reduced excitation amph- 
tudes r^'^{J) and rj'"-"'{Jab,J)- In this coupled scheme, 



the eigenvalue equation in Eq. ([T6| is solved separately 
for each set of quantum numbers { J'^, T^}. 

A simple estimate shows that a model space of Uq + n^ 
single-particle states consists of only (no+n„)^/'^ j-sheUs. 
Thus, the entire computational effort is approximately 
reduced by a power 2/3 within the spherical scheme com- 
pared to the TO-scheme. We have derived and numerically 
implemented the spherical scheme within the CCSD ap- 
proximation. We tested that our TTi-scheme code and 
the spherical code give identical results for several test 
cases. The spherical code permits us to reach much 
larger model spaces, and we are able to achieve satisfac- 
tory convergence even for "hard" interactions such as the 
"bare" N^^LO interaction. In Ref. [8| we used the spher- 
ical coupled-cluster code to calculate the ground states 
of medium-mass nuclei like ^°'^®Ca and "^^Ni within the 
CCSD approximation and starting from bare chiral in- 
teractions. In a model space of 15 major oscillator shells, 
the results were reasonably well converged. Recently, 
we have also implemented the A-CCSD(T) approach in 
an angular-momentum-coupled scheme. In Ref. [23|, we 
calculated the ground states of the oxygen isotopes us- 
ing chiral interactions in model spaces comprising up 
to 20 major oscillator shells within the A-CCSD(T) ap- 
proach. In such large spaces we were even able to con- 
verge the ground state energies for a chiral interaction 
with a 600 MeVc~^ momentum cutoff. Likewise, the in- 
clusion of continuum scattering states - necessary for a 
description of halo states in weakly bound nuclei - yields 
large model spaces that can be treated within he spherical 
scheme [l^l- Details of the angular- momentum-coupled 
A-CCSD(T) approach are given in the Appendix. 



Here E{13) is the ground-state energy of the Hamilto- 
nian H + /SB. Unfortunately, the relation (P^ does not 
hold exactly within the coupled-cluster approach since 
this method is not variational when the cluster operator 
is truncated. We will nevertheless base the computa- 
tion of expectation values on Eq. ((23|). Experience shows 
that this approach is approximately valid if the relevant 
particle-hole excitations are incorporated. In practice, 
the differential quotient (j23p is numerically implemented 
as a difference quotient with /3 « 0.01. This requires us to 
perform two calculations (one for zero and one for small 
nonzero /?) . For relatively "hard" interactions such as the 
"bare" N^LO, we also perform two Hartree-Fock calcu- 
lations and thus employ slightly different single-particle 
bases in these two cases. While such a procedure is un- 
necessary for methods that fulfill the Hellmann-Feynman 
theorem, it is important within our approach. We note 
that this approach most closely reflects the physical sit- 
uation of an expectation value measuring the system's 
response to a perturbation 



III. TREATMENT OF THE CENTER-OF-MASS 
PROBLEM 

In this section, we demonstrate that the coupled- 
cluster wave function factorizes to a very good approxi- 
mation into a product of a center-of-mass wave function 
and an intrinsic wave function. We will present a simple 
procedure that checks this factorization and give an es- 
timate for the degree of the achieved factorization. This 
section significantly expands on the short demonstration 
of the factorization in Ref. ^] . 



D. Expectation values in coupled-cluster theory 

For expectation values and the computation of ground 
state properties other than the energy, we utilize the right 
and left eigenvectors of the similarity-transformed Hamil- 
tonian H and compute the one- and two-body reduced 
density matrices. 



Ppq = {<Po\^e alaqe |(/)o) 



Ppqrs 



ho\Ae~-^ alalarUse^ |(/)o) 



(21) 
(22) 



Expectation values can then be computed by expressing 
the operator of interest in terms of the density matri- 
ces. The Hellmann-Feynman theorem is another route 
to the computation of expectation values. In this case, 
we consider a response of the ground-state energy to a 
perturbation caused by a given operator of interest. The 
Hellmann-Feynman theorem expresses the expectation 
values of an observable B as 



(B) = 



dE{(3) 



dp 



= 



(23) 



A. Statement of the center-of-mass problem 



Let us consider the nuclear ^-body Hamiltonian 



ha^y: 



i=i 



2m 



j<k 



(24) 



Here, 1^ is a two-body operator that is invariant under 
rotations and translations. Thus, the total momentum 
and angular momentum are conserved quantities. It is 
advantageous to separate the Hamiltonian into an intrin- 
sic Hamiltonian Hin and the center-of-mass Hamiltonian 
(i.e., the kinetic energy Tcm of the center of mass) as 



H — Tcm 

= T, 



A 



m 



cm I 7 , 

]<k 



2mA 



V{j,k) 



(25) 



Note that the intrinsic Hamiltonian Hm does not depend 
on the center-of-mass coordinate. Thus, we could add an 
arbitrary operator Hem of the center-of-mass coordinate 
to the intrinsic Hamiltonian Hm without changing the 



intrinsic properties of the resulting Hamiltonian. In other 
words, the eigenfunctions ipA of the A-body Hamiltonian 
Hin+Hcm are products of an intrinsic eigenfunction ipi^ of 
the intrinsic Hamiltonian Him and a center-of-mass wave 
function -tpcm that is the eigenfunction of the center-of- 
mass Hamiltonian iJcm, that is, 



i'A = V'cmV'ii 



(26) 



The corresponding energy is the sum Ea — -Ecm + E\^\ of 
the center-of-mass energy and the intrinsic energy. The 
question thus arises which operator i?cm to choose. Let 
us consider the NCSM Q as an example. Here, one works 
in a complete Nhuj space consisting of a basis of all A- 
particle Slater determinants of oscillator states with fre- 
quency u) and total excitation energy not exceeding Nhjj. 
One can add the center-of-mass Hamiltonian 



Hcm{^) — Tc, 



-mAuj'^R, 
2 



2 
cm 



— tub 

2 



(27) 



with cj = cj to the intrinsic Hamiltonian _ffin, and ob- 
tains a factorizcd ground-state wave function where the 
center-of-mass wave function i/jcm is a Gaussian with fre- 
quency w. The (truncated) coupled-cluster method is un- 
able to employ an Nhiu space, and it is arguably the best 
idea to completely remove any reference to the center- 
of-mass coordinate. Thus, we solve an intrinsic Hamil- 
tonian TJjn that depends on ^ — 1 independent coordi- 
nates in a Hilbert space spanned by wave functions of 
A coordinates. Two comments are in order. First, the 
truncated coupled-cluster method is not an exact solu- 
tion of the A-body problem, but rather a very efficient 
approximation. Thus, it is not guaranteed a priori that 
the coupled-cluster wave function exhibits the factoriza- 
tion (P5|. Recall that any A-body wave function can be 
expanded as 



V'^ = Vs,v^M^ 



i>i 



(28) 



where X^i ■^f ~ ^ from normalization, and we assume 



that the non-negative weights Sj are ordered in decreas- 
ing order. Only if all but one of the weights Sj van- 



ishes, does the factorization \2Q\ take place; otherwise, 
the factorization might be mildly or strongly violated, 
depending on the size of the weights Sj. Second, if a fac- 
torization of the coupled-cluster wave function takes (ap- 
proximately) place, the form of the center-of-mass wave 
function -i/jcm corresponding to the largest weight si has 
to be determined. In what follows we will see that the A- 
body coupled-cluster wave function factorizes to a good 
approximation (i.e., si ~ 1), and that the center-of-mass 
wave function is the Gaussian ground state of the center- 
of-mass Hamiltonian (l?71) for a yet-to-be determined fre- 
quency U). 



B. Approximate separation of the center-of-mass 
wave function in coupled-cluster calculations 

We consider the nucleus ^^O and employ the low- 
momentum interaction T^ow fc with a smooth momen- 
tum cutoff A = 1.8 fm~^ derived from Machleidt's and 
Entem's chiral N'^LO interaction. Figure [T] shows the 
ground-state energy as a function of the oscillator spac- 
ing hhj of the underlying oscillator basis in a model space 
of nine oscillator shells. These results are obtained from a 
CCSD calculation within a spherical Hartree-Fock basis. 
The energy is well converged as it depends very weakly 
on the model-space parameters. We also computed the 
ground-state expectation value 



Er 



{H{u)) 



(29) 



of the center-of-mass Hamiltonian ((27)) for a frequency 
w = w and show the result in the inset of Fig. [TJ This 
expectation value is generally not zero. This indicates 
that the coupled-cluster wave function is in general not 
an eigenstate of the center-of-mass Hamiltonian H(uj). 
However, there seems to be little correlation between 
the ground-state energy and the expectation value E(ll)), 
and the latter vanishes approximately in a model space 
with hoj Ri 20 MeV. Thus, at this frequency, the coupled- 
cluster wave function is approximately the ground state 
of the center-of-mass Hamiltonian H{uj). As a check, 
we fix /iw = 20 MeV, consider the Hamiltonian H = 
-ffin + /3H{uj), and compute its ground-state energy as 
a function of the parameter /3. The result is shown in 
Fig. [21 Clearly, the ground-state energy is rather insen- 
sitive to /3 and varies by only 15 keV as /3 is increased 
from zero to one. 
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FIG. 1: (Color online) Ground-state energy (within CCSD) 
of ^"O with a low-momentum interaction as a function of the 
oscillator spacing hoj. The model space consists of nine ma- 
jor oscillator shells. Inset: Expectation value _Ecm(^) of the 
center-of mass Hamiltonian with the standard frequency de- 
pendence. 
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FIG. 2: (Color online) Ground-state energy (within CCSD) 
of ^®0 using a low- momentum interaction with the center- 
of-mass Hamiltonian /3-ffcm added. Calculations employ nine 
major oscillator shells and hu — 20 MeV. As /3 is varied be- 
tween and 1, the ground state energy changes by about 
15 keV. 



The factorization of Eq. ((26|) is thus achieved in a 
model space with fku ^ 20 MeV. Let us assume that 
such a factorization takes place generally, and that the 
corresponding center-of-mass wave function is a Gaus- 
sian with unknown frequency w that might differ from 
the frequency of the underlying oscillator basis. Thus, 
we assume that the coupled-cluster wave function is the 
ground state of the center-of-mass Hamiltonian ([27]) for 
a suitable frequency uj. To determine this frequency we 
employ the identity 



3 cj^ 

Hc-m{uj) + -hoj- Tern = "7^ ( i?cm(t^) 



— huj — Tr, 

2 



and take its expectation value. We seek (and insert) 
-E'cm(w) = 0, employ the relation (Tcm) — |to valid for 
Gaussians, insert the already computed expectation val- 
ues Ecmi^), and solve for the unknown frequency w. This 
yields the two possible solutions. 



hlJ = fuv+ -£'cm(w) 



±^^iE,miuW + ^hluEcra{L0) . (30) 

We compute the ground-state expectation values Ecmi^jj) 
for these two frequencies and find that one expectation 
value is typically very close to zero. Figure [3] shows 
that the small expectation value essentially vanishes for a 
large range of frequencies of the underlying oscillator ba- 
sis in a model space of 13 major oscillator shells. Closer 
inspection shows that the expectation value is about 
-Ecm(w) « — 10 keV. Recall that coupled-cluster theory 
is non-variational (as the similarity-transformed Hamil- 
tonian is non-Hermitian) , and such a small negative ex- 
pectation value is certainly tolerable for the non-negative 
operator (P7|) . Figure |3] also shows that the frequency of 
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FIG. 3: (Color online) ^^0 Ground-state expectation value 
(within CCSD) of the generalized center-of-mass Hamiltonian 
Hci-n{ij^) and of the kinetic energy Tcm as a function of the 
oscillator spacing huj. The model space consists of thirteen 
major oscillator shells. 



the Gaussian center-of-mass wave function stays approx- 
imately constant fiw w 20 MeV for a large range of fre- 
quencies w of the underlying oscillator basis. As a final 
check, we also computed the ground-state expectation 
value of the center-of-mass kinetic energy. For a Gaus- 
sian, this expectation value fulfills (Tcm) = f^^- Figure^] 
shows that this relation is indeed obeyed. 

Let us turn to the "bare" N^LO interaction [Hill]. 
This interaction has a substantially higher momentum 
cutoff than the previously employed low-momentum in- 
teraction. For converged results, we need to employ 
larger model spaces and also triples clusters. We per- 
formed coupled-cluster calculations of the ground states 
of ^^O and "^He. Figure S] shows the ground-state en- 
ergy of ^®0 (lower panel), the ground-state expectation 
value £'cm(w) (upper panel) of the center-of-mass Hamil- 
tonian (1271) . and the corresponding frequency uj (middle 
panel) as a function of the frequency uj of the underlying 
oscillator basis. One sees that the ground-state energy 
is well converged and displays only a weak dependence 
on UJ. Similar comments apply to the frequency w. The 
expectation value Ecmiuj) is small, i.e., |ii'cm('^)| ^ ^^ 
but nonzero. For smaller values of the frequency w of the 
underlying oscillator basis, we are at the limit of well- 
converged results, and Ecm even becomes negative. This 
shows that the coupled-cluster ground state is not an 
exact eigenstate of the center-of-mass Hamiltonian (|27p . 
Let us also assume that the factorization is not perfect 
(though this is not implied by a nonzero value of Ecm)- 
We want to estimate the level of admixture of center- 
of-mass excitations. A simple and conservative estimate 
based on perturbation theory shows that the admixture 
is essentially Ecm/fiui « 5%. We also computed the in- 
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FIG. 4: (Color online) Bottom: Ground-state energy of ^"O 
within the A-CCSD(T) approximation as a function of the 
frequency hw of the underlying oscillator basis. Middle: Re- 
lation between the frequency uj and the frequency uj of the 
underlying oscillator basis. Top: Expectation value _Bcm('i) 
of the center-of-mass vs. hw. 



trinsic point radius r from 



(31) 



l<j<k<A 



and show the result in Fig. [5] This calculation is again 
based on the triples-corrected A-CCSD(T) approxima- 
tion, and it still exhibits a weak model-space dependence. 
It provides us with a further test of the wave-function 
factorization, since we can employ the center-of-mass co- 
ordinate i?cm and rewrite Eq. (PT|) as 



1 ^ 



1 ^ 

((f,-i?cmf) = -5](rf)-(i?L) 



(32) 



j=i 



For a Gaussian corresponding to the frequency uj, we have 
- 3fc,r, 



^mACj'^{R^^ — |fta), and the intrinsic point radius can 



thus be computed from the expectation value of a one- 
body operator. The result from this calculation is in very 
good agreement with the result obtained from Eq. ([3T|). 
as shown in Fig. [5l This suggests that the factorization 
of the wave function might be better than expected from 
the calculation of the ground-state expectation value ([^^ 
depicted in Fig. |4l 

At the moment, we have no analytical insights into 
the observed factorization and the Gaussian shape of the 
center-of-mass wave function. We believe that the factor- 
ization itself is not entirely surprising. While our finite 
basis is not complete in a mathematical sense, it is suffi- 
ciently complete to describe low-energy nuclear structure, 
i.e., it contains momentum modes that exceed the cutoff 
of the interaction, and it is sufficiently extended in posi- 
tion space to accommodate a quantum object of the size 
of the nucleus. Thus, one expects the basis to capture 
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FIG. 5: (Color online) RMS radii of ^"0 within the A- 
CCSD(T) approximation using the intrinsic point radius and 
subtracting (-Rem) assuming a Gaussian for the center of mass. 



the relevant physics, and a well-factorized ground state is 
not too surprising, even if the employed many-body basis 
functions do not individually reflect this factorization. 



C. Center-of-mass problem for open-shell nuclei 

Odd-mass nuclei that differ from a nucleus with closed 
subshells by the addition or removal of one nucleon can 
also be computed efficiently within the coupled-cluster 
method. For the calculation of ^^O, we employ the 
particle-attached equation-of- motion method (J16p . which 
is based on describing ^^O as a nucleon (i.e a superposi- 
tion of Ip and 2p-lh excitations) upon the ground state 
of ^^O. We first solve the coupled-cluster equations (|H) for 
the "mass-shifted" nucleus ^^O, i.e., for the computation 
of the ^^O ground state, we employ the Hamiltonian ([M)) 
with mass number A — 17. In a second step, we describe 
the ground and excited states in ^^O in terms of the ex- 
citations (fT5|) upon this ^^O ground state. For an SRG 
interaction with cutoff A = 2.8 fm~^, the energies of a 
few low- lying states are shown in Fig. [S] as a function of 
the oscillator frequency in a model space with TV = 12. 
These states are single-particle states and can be well 
computed within our approach. Figure [7] shows the en- 
ergy expectation value i?cm(w) (defined in Eq. (f^ ) for 
the three low-lying single-particle states as a function of 
the oscillator frequency lo of the underlying model space. 
These expectation values are small for a; sa 16 MeV, but 
large for other parameters. This shows that the center- 
of-mass wave function is generally not a Gaussian with 
frequency lo. No further conclusion can be reached from 
this expectation value. In particular, this is no evidence 
that the center-of-mass wave function does not separate. 

Let us again assume that the coupled-cluster wave 
function factorizes into a Gaussian (with frequency uj) 
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FIG. 6: (Color online) Energies of low-lying J'^ = 1/2+ , J'' = 
3/2+, J^ = 5/2+ states in ^^O based on an SRG interaction 
with cutoff A = 2.8 fm"^ in a model space consisting of Af + 
1 = 13 oscillator shells versus the oscillator frequency. 
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FIG. 8: (Color online) Frequency Cj of the (approximately) 
Gaussian center-of-mass wave function for the low-lying J^ = 
1/2+ , J'' = 3/2+, ,r = 5/2+ states in "O versus the oscilla- 
tor frequency. 
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FIG. 7: (Color online) Expectation values £'cm(i^) of the 
center-of-mass Hamiltonian 7Jcm(w) (see Eq. (|27|l ) for the low- 
lying ,r = 1/2+ , .r = 3/2+, r = 5/2+ states in ^^O versus 
the oscillator frequency. 
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FIG. 9: (Color online) Expectation values _Ecm(i^) of the 
center-of-mass Hamiltonian Hcmi^) (see Eq. (|27p ') for the low- 
lying J"" = 1/2+ , J"" = 3/2+, J'" = 5/2+ states in ^'^O versus 
the oscillator frequency. Here, the frequency di = d;(cj) em- 
ployed in the center-of-mass Hamiltonian is as depicted in 
Fig. [31 We have Ecm{ijj) <g hw and thus a very good separa- 
tion of center-of-mass excitations. 



for the center of mass and an intrinsic "wave function, 
and let us determine Co from Eq. (pO]). Figure |8]sho-ws the 
resulting frequency for the three states we computed in 
^^O. The corresponding expectation values £'cm('^) de- 
fined in Eq. ([2^ are sho-wn in Fig. [HI as a function of 
the oscillator frequency lo of the underlying model space. 
"We employ uj = Cj{uj) as depicted in Fig. [H As evident 
from Figs. |8] and El -we have £'cm(w) <C hCo. Thus, the 
computed states in ^^O also exhibit a Gaussian center-of- 
mass "wave function to a very good approximation (ad- 
mixtures of higher hw oscillator excitations are of the 
order of Ec^{6j) / huj « 1-2%). This implies that the fac- 
torization is at least of the same degree of quality. 



D. Center-of-mass problem in a toy model 

To shed some Ught onto the observed factorization, we 
consider a simple problem of t'wo interacting particles in 
one spatial dimension. We choose a two-body interaction 
of the form V{x) = — Vb exp(— (a;//)^), where x = {xi — 
X2)/V2 is (up to a factor v^) the relative coordinate of 
the two-particle system, and I is a length scale. "We thus 
consider the intrinsic Hamiltonian 
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P 
2m 



Vbexp(-(x/0') 



(33) 
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with m being the mass and p the relative momentum 
p = {pi — P2)/V^. We choose a basis consisting of prod- 
ucts ^rn{xi/l)^n{x2/l) of oscillator wave functions $fc, 
and choose < m,n < N. Thus, the basis is not a com- 
plete Nhu space, and we solve an intrinsic Hamiltonian 
that depends on the relative coordinate in a model space 
consisting of single-particle coordinates. We consider dif- 
ferent interaction parameters Vq, and set the oscillator 
length of our basis equal to the scale I of the Gaussian 
interaction. Note that the resulting ground state has an 
extension that, depending on the strength Vq of the inter- 
action, differs considerably from I. Indeed, approximat- 
ing the interaction by a parabola at its minimum shows 
that the corresponding frequency is il/uj — y^2Vo/{fiuj). 
Again, we find that the ground-state wave function has 
a Gaussian shape in the direction of the center-of-mass 
coordinate {xi +X2)/V2. Figure [TOl quantifies this state- 
ment. The circles show the relative error of the ground- 
state energy (obtained from comparing the result in a 
model space of A^ = 8 oscillator shells with the result in 
A^ = 16 oscillator shells) as the strength Vq of the interac- 
tion is varied. The squares (diamonds) show to what ex- 
tent the ground-state wave function factorizes in a model 
space consisting of A^ = 8 (A^ — 16) oscillator shells. 
The residual 1 — s^ is clearly very small, and it decreases 
with increasing size of the model space. The triangles 
show the ground-state expectation value of the center-of- 
mass Hamiltonian (|27p , normalized by the corresponding 
spacing huj. Thus, the center-of-mass wave function fac- 
torizes, and it is to a high accuracy a Gaussian. The 
result obtained for this simple Hamiltonian lends further 
support to the results obtained in the coupled-cluster ap- 
proach. It thus seems that one might obtain a very accu- 
rate factorization in sufficiently large model spaces even 
without employing an Nhuj space. Again, this factoriza- 
tion is not surprising since the model space becomes more 
and more complete as it increases in size. The emergence 
of a Gaussian center-of-mass wave function, however, is 
remarkable and not yet understood. 

The center-of-mass problem in coupled-cluster calcu- 
lations has also recently been addressed by Roth and 
coworkers 48]. These authors consider a Hamiltonian 
of the form H ~ Hin + (3Hcm{^), that is, the center-of- 
mass Hamiltonian (l?fl) with a frequency fixed to that of 
the employed harmonic oscillator basis is added to the in- 
trinsic Hamiltonian. This method, pioneered by Gloeck- 
ner and Lawson [49|, aims at shifting spurious states up 
in the spectrum. Roth et al. found that the coupled- 
cluster wave function for such a Hamiltonian does not 
factorize into a Gaussian with frequency uj for center- 
of-mass wave function and an intrinsic wave function. 
Two comments are in order. First, we believe that this 
approach misses the point as it requires perhaps more 
than the coupled-cluster approximation needs to deliver. 
For atomic nuclei, there is no center-of-mass interaction, 
and the requirement that the wave function factorizes 
in presence of a center-of-mass interaction is more than 
one needs. Some time ago, McGrory and Wildenthal also 
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FIG. 10: (Color online)Results for the toy model (|33|l as a 
function of the strength of the interaction Vq. Circles: Esti- 
mate of relative error A_E/_E of the ground-state energy by 
comparing a calculation in A = 8 oscillator shells with a cal- 
culation in A = 16 shells. Squares and diamonds: Fraction 
1 — Si of the ground state wave function that is not factor- 
ized in an intrinsic state and a center-of-mass state, obtained 
from the singular value decomposition (|28|) in a model space 
of A'^ = 8 and A'^ = 16 oscillator shells, respectively. Triangles: 
The ground-state expectation value Ec-m{ijj) of the harmonic 
oscillator center-of-mass Hamiltonian is much smaller than 
the energy hw of spurious center-of-mass excitations, indicat- 
ing that the center-of-mass wave function is approximately a 
Gaussian. 



pointed out that the method by Gloeckner and Lawson is 
not fully appropriate in truncated model spaces that are 
not complete Nfiw spaces [50l| . For nuclei, one only needs 
that the wave function of the intrinsic Hamiltonian, when 
computed in a Hilbert space of A coordinates, exhibits 
a factorization. Second, we caution that the conclusions 
by Roth and coworkers are based on results that are not 
fully converged with respect to the size of the model space 
(see Sect. Em below). 



The results of this section suggest the following proce- 
dure when dealing with model spaces that are not com- 
plete Nhw spaces: (i) Compute the spectrum of the in- 
trinsic Hamiltonian H\^ in as large a model space as con- 
veniently possible, (ii) Check a posteriori whether the 
resulting center-of-mass wave function is a Gaussian by 
computing the expectation value Ecia., and determine 
the corresponding frequency Q. (ni) The ratio £^cm/(fi^) 
serves as a conservative estimate for the quality of the 
factorization. (Strictly speaking, this ratio measures to 
what extent the center-of-mass wave function deviates 
from a Gaussian with frequency uj.) 



11 



IV. COUPLED-CLUSTER RESULTS FOR 

MEDIUM-MASS NUCLEI WITH CHIRAL 

INTERACTIONS 

In this section, we compute binding energies and radii 
of several doubly magic nuclei within the CCSD and the 
A-CCSD(T) approximation, and employ chiral NN inter- 
actions. 

Our single-particle wave functions are eigenfunctions 
of the harmonic oscillator and characterized by the fre- 
quency ui. Our model space consists of spherical os- 
cillator wave functions with radial quantum number n 
and angular momentum /, and we include single-particle 
states with 2n + l < N and I < 10. Fully converged ener- 
gies have to be independent of the parameters N and lu 
of our single-particle basis. In practice, we cannot go to 
infinitely large spaces, but this is also not necessary for 
the description of low-energy properties of finite nuclei. 
Our basis needs to be complete in the following sense. 
It must be sufficiently extended in momentum space to 
resolve the cutoff A^ of the employed interaction, and its 
extension in position space must be such that a nucleus 
of radius R literally fits into the basis. Let us estimate 
the parameters required for a model space of oscillator 
functions. 

An oscillator basis consisting of N shells at a frequency 
uj resolves the high momentum cutoff A if the inequality 



2\2 



Nhuj > 



h^x 



(34) 



is fulfilled. Likewise, the basis has to be sufficiently ex- 
tended in position space to describe a nucleus with a 
radius R. For this, the inequality 



huj <N- 



iR^ 



(35) 



needs to be fulfilled. In other words, an oscillator basis 
with oscillator length /ho = \/h/ (muj) exhibits the in- 
frared cutoff l^l I^N and the ultraviolet cutoff I^ioVN- 
Thus, converged results require model space parameters 
(iV, w) that fulfill the inequalities (l34l) and ((35)) . and the 
results will then be insensitive to the specific values of 
the parameters. Note that the simultaneous fulfillment 
of Eq. dnS]) and Eq. dH]) requires N > XR. (Our calcu- 
lations presented in this paper show that these approxi- 
mate relations are reasonable estimates.) Note also that 



hw 



mR^ 



(36) 



for a minimum model space with N « XR. Thus, the 
well-known estimate huj « 42/^^/'^ MeV is only valid for 
small cutoffs X Ri kp that are close to the Fermi mo- 
mentum kp. These considerations show that much is to 
be gained from low-momentum interactions, and - con- 
versely - that "bare" chiral interactions with A = A^ 
require very large model spaces. 



Let us consider the oscillator basis as an example. A 
model space of N oscillator states contains about il ~ 
N^/3 single-particle states. We obtain converged results 
for ^^O and the interactions from chiral EFT in TV « 15 
shells. In this model space there are ni/{{il — Z)\Z\) k, 
10^° Slater determinants for the protons alone. Thus, the 
resulting model space for protons and neutrons is far out 
of reach from diagonalization methods. 

We employ the chiral NN interaction at order N'^LO, 
i.e., the NN interaction is included up to the order 
(Q/A^)^, while 3NF at this order and interactions at 
higher order are neglected. We perform a spherical 
Hartree-Fock calculation and transform the Hamiltonian 
to the Hartree-Fock basis. The CCSD equations are 
then solved in the spherical Hartree-Fock basis, and the 
Hartree-Fock state is taken as the reference state for the 
coupled-cluster method. Note that the chiral N^LO inter- 
action is quite "hard" . Within the Hartree-Fock approx- 
imation one does not even obtain bound nuclei. Nev- 
ertheless, the solution of the non-linear coupled-cluster 
equations yields rather well-bound nuclei. 

Fig. [Tl] shows the ground-state energy of '*He in a 
model space of A^ = 18 (19 major oscillator shells). In 
such a large model space, the results are virtually inde- 
pendent of the frequency uj over a wide range. The CCSD 
result for ^He deviates from the virtually exact result 
from the solution of the Faddeev-Yakubowsky (FY) equa- 
tions by about 6%. This difference is due to the omission 
of three- and four-particle clusters. The more accurate 
A-CCSD(T) method includes three-body clusters approx- 
imately and overshoots the (FY) result by about 1.5%. 
This ovcrbinding is due to the non- variational character 
of the coupled-cluster method and is not really a concern 
due to the accuracy of the method. The experimental 
ground-state energy of ^He \s E = — 28.3MeV, and the 
additional binding must be attributed to the 3NFs at 
order N'^LO; other high-order terms in the chiral Lan- 
grangian play only a very small role in He. 

We also computed the radius of the alpha particle 
within the Hellmann-Feynman approach. The results 
shown in Fig. 1121 exhibit a very weak dependence on the 
oscillator frequency, particularly within the A-CCSD(T) 
approximation. Within this approach, the radius de- 
viates by about 0.01 fm from the NCSM results Hj], 
and about 0.015 fm from the virtually exact result from 
the hyperspherical harmonics method employed by the 
Pisa group [4J . We believe that the results presented for 
the alpha particle manifest the high degree of accuracy 
that the coupled-cluster approximation exhibits. Note 
that the alpha particle also exhibits the approximate fac- 
torization into a Gaussian center-of-mass wave function 
and an intrinsic wave function. Here, the frequency of 
the Gaussian is to « 19 MeV while the expectation 
value of the center-of-mass Hamiltonian is much smaller 
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0.3 MeV. 



We next turn to the calcium isotopes ''"Ca and ^®Ca. 
Figureinishows the CCSD and A-CCSD(T) ground state 
energy for *°Ca as a function of the oscillator frequency 
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FIG. 11: (Color online) CCSD and A-CCSD(T) ground state 
energy of *He using a chiral NN interaction at order N'^LO 
in 19 major oscillator shells as a function of the oscillator 
spacing hu, compared to virtually exact results from Faddeev- 
Yakubowsky calculations. 



FIG. 13: (Color online) A-CCSD(T) and CCSD results for 
^"Ca from the chiral NN interaction at order N"^LO as a func- 
tion of the oscillator spacing &j and the size of the model 
space. 
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FIG. 12: (Color online) *He RMS radii calculated using the 
Hellmann-Feynman theorem within the CCSD and the A- 
CCSD(T) approximations using a chiral NN interaction at 
order N"^LO with cutoff 500 MeVc~^. Calculations employed 
19 major oscillator shells, and the results are plotted as a 
function of the oscillator spacing huj. 



hw and the size of the model space TV (here the num- 
ber of major shells is A^ -f 1). For hw = 32 MeV and 
N + 1 = 15 sheUs the A-CCSD(T) ground state energy 
is —345.074 MeV. We also performed a calculation in 
iV+1 ^ 19 shells at ^iw = 32 MeV, yielding a A-CCSD(T) 
ground state energy of —345.781 MeV. This shows that 
our results are converged within 1-2 MeV with respect 
to the size of the model space. Within this approxima- 
tion, one overbinds ''"Ca by about 3 MeV. Compared to 
the previously published CCSD results, the triples cor- 
rections add more than 30 MeV of binding. 

Figure [H shows the A-CCSD(T) resuhs for ^SCa. The 



results are converged to within about 2 MeV in a model 
space of A'' -|- 1 = 19 shells, and only exhibit a weak 
model-space dependence in the largest model space. The 
approximation with triples yields about an additional 
40 MeV of binding when compared to the CCSD results. 
Let us assume that the difference between our calcula- 
tions and the experiment is mainly due to the omitted 
3NFs. This suggests that the three-nucleon forces (3NF) 
will exhibit an interesting isospin dependence. For *'^Ca, 
only small (repulsive) contributions from the 3NF are ex- 
pected, while about 20 MeV of attraction is needed for 
^®Ca. A similar picture arises when comparing ^^O with 
^^O. For the former, about 0.41 MeV per nucleon in bind- 
ing are missing, while the latter lacks about 0.82 MeV per 
nucleon in binding [27[. 
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FIG. 14: (Color online) A-CCSD(T) and CCSD results for 
^*Ca from the chiral NN interaction at order N^LO as a func- 
tion of the oscillator spacing huj and the size of the model 
space. 
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At this point, it is interesting, to compare the ener- 
gies from the Hartree-Fock calculation, the CCSD calcu- 
lation, and the A-CCSD(T) approximation. This com- 
parison will permit us to estimate the precision of the 
coupled-cluster method. Table U shows that CCSD ac- 
counts for about 90% of the correlation energy, while the 
A-CCSD(T) approximation yields about 10%. This hier- 
archy has been observed in quantum chemistry as well, 
see for example Ref . [3j| . For the alpha particle we know 
that four-particle clusters account again for about 10% 
of the triples correction. Thus, it makes no sense to arti- 
ficially increase the precision of our results by turning to 
even larger model spaces. The error estimates due to the 
finite model space are of the same order as the estimates 
due to omitted four-body clusters. Note that the results 
in Table |T] are also consistent with the size extensivity of 
the employed coupled-cluster methods, as the deviation 
AE from data is approximately linear in mass number, 
i.e., AE/A is approximately constant over a considerable 
range. 



(1^/2, and (1^/2 proton and neutron single-particle states in 
the closed-shell oxygen isotopes ^^O, ^^O, '^^O, and ^*0. 
Here, we define the single-particle energy of the state J'^ 
as the difference in binding energy between the A ± 1 
nucleus and the closed-shell nucleus A, i.e., 



A±lf TW\ 



E,pir) = E^^'ij 



E' 



(37) 
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The evolution of nuclear shell structure for neutron- 
rich isotopes is of great theoretical and experimental in- 
terest, see for example the recent review [54]. Some of 
the traditional magic numbers of the nuclear shell model 
might fade away as one moves away from the valley 
of beta stability, while new magic numbers emerge in 
neutron-rich nuclei. The evolution of shell structure in 
the isotopes of ox yge n has recently received considerable 

ferimental [Sa - ISSf and theoretical attention [23, [59l - 
Considerable shell gaps have been observed in ^^O 
and ^^O [5a,|57|, leading to the interpretation of new 
magic numbers at (Z = 8, iV = 14) and (Z = 8, iV = 16), 
respectively. Hoffman et al. [58] found that ■^''^O is a res- 
onance and unstable towards one-neutron emission. The 
instability of ^^O is closely related to the location of the 
1^3/2 single-particle shell, and it is clear that the evolu- 
tion of the d^/2 shell will decide whether ^^O or ^^O is 
the most neutron-rich stable isotope of oxygen. 

The underlying microscopic mechanisms for the evolu- 
tion of nuclear shell structure is not yet entirely under- 
stood. One might expect that both three-nucleon forces 
and coupling with the scattering continuum will play sig- 
nificant roles in predicting the limits of nuclear stability 
and the shell evolution towards the drip line. Zuker [GJ] 
suggested that three-nucleon forces modify the monopole 
terms of microscopically derived shell-model interactions. 
Indeed, Otsuka et al. [62i] found within the sd shell model 
that 3NFs will add repulsion between the (^5/2 and ^3/2 
orbitals, thus making ^"^O the heaviest bound isotope 
of oxygen. Coupled-cluster calculations of the isotopes 
i6,22,24,28q j-gyeal that three-nucleon forces (or more com- 
plicated many-body forces) play a non-negligible role in 
the determination of the drip-line for oxygen isotopes. 

In this section, we study the evolution of shell struc- 
ture in the isotopes of oxygen and fluorine based on chiral 
NN interactions [13|] . We will omit continuum effects and 
3NFs from our study. Recall that chiral 3NFs consist of 
a long-ranged two-pion exchange, a mid-range one-pion 
exchange, and a contact term [GSj . Variation of the cutoff 
will induce short-range 3NFs that are in their structure 
identical to the chiral 3NF [gi, [13 ■ Thus, we will be 
able to probe the effect of short-ranged 3NFs a posteri- 
ori by cutoff variation. Continuum effects are, of course, 
expected to be important for nuclei close to the drip 
line 18, 68, 69]. However, it will turn out that most of 
the computed single-particle states are fairly well bound 
(likely due to the omission of 3NFs). In this situation, 
the continuum plays a smaller role. 

Figures [TS] and [11] show the effective single-particle 
in oxygen and fluorine are usually interpreted as the si/2, energies for the neutron and proton Sl/2^ ^3/21 and ^5/2 



TABLE I: Cround-state energies per nucleon E/A and devi- 
ation AE/A from experiment for doubly magic nuclei within 
the CCSD and A-CCSD(T) approximations. All energies are 
in units of MeV. 

It is also very interesting to compare our results with 
the results [52| by Fujii et al. obtained within the Uni- 
tary Model Operator Approach [5^. For ^^O, Fujii et 
al. report binding energies of 6.62 MeV per nucleon and 
7.47 MeV per nucleon employing two-body clusters and 
three-body clusters, respectively. The close agreement 
between our results presented in Table U and those of 
Ref. [52| demonstrate that different ah initio methods 
are setting reliable benchmarks for increasingly heavier 
nuclei. 



V. SHELL EVOLUTION WITH CHIRAL NN 
INTERACTIONS 

In this section we investigate the evolution of the low- 
lying positive parity states 1/2 , 3/2 , and 5/2 in 
the oxygen and fluorine isotope chains using the chiral 
nucleon-nucleon interaction at N^LO. Our calculations 
do not include three-body forces, but in order to probe 
the effects of omitted short-ranged many-body forces we 
renormalize the "bare" chiral interaction using the simi- 
larity renormalization group method and study the reso- 
lution scale dependence on the calculated energies. The 
low-lying states J'^ = 1/2+ , J^ = 3/2+, and r = 5/2+ 
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states in ^^O, ^^O, ^''O, and ^^O. The calculations em- 
ployed a model space consisting of A^ + 1 = 15 oscilla- 
tor shells at a fixed oscillator frequency hoj — 28 MeV. 
The proton si^, (1^/21 E'-nd the d^/2 single-particle states 
were computed within the PA-EOM-CCSD approach. 
The neutron single-particle states resulted from either 
the PR-EOM-CCSD or the PA-EOM-CCSD method (See 
subsection III Bl for more details), depending on whether 
the considered single-particle state is a hole or a particle 
state of an isotope of oxygen with a closed subshell. 




FIG. 15: (Color online) Single-particle energies (Eq. (|37|) ) of 
the Si/2, ^3/2, and the dg/2 single-particle states in the oxygen 
isotopes using the "bare" chiral NN interaction. The calcu- 
lations employ A'' = 14 shells for fixed hw = 28 MeV. The 
^5/2 and Si/2 experimental single-particle energies are shown 
as dashed lines. 
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FIG. 16: (Color online) Single-particle energies (Eq. (|37p l of 
the Si/2, 0(3/25 and the dc,p single-particle states in the flu- 
orine isotopes using the "bare" chiral NN interaction. The 
calculations employ Ai" = 14 shells for fixed fej = 28 MeV. 
The experimental ^5/2 single-particle energiy is shown as a 
dashed line. 

Several interesting features can be extracted from the 
results. The computed neutron single-particle energies 
(shown in Fig. I15p do not reproduce the experimentally 
observed shell gaps in ^^O and ^"'O. Furthermore, the 



chiral NN interactions incorrectly yield a bound ^^O j58|. 
This behavior has also been seen in shell-model calcula- 
tions of the oxygen isotopes using microscopically derived 
shell- model interactions (see for example Ref. [62'|). For 
the fluorine isotopes (shown in Fig. [T6|) . the 1^5/2 state 
follows the experimental trend but lacks binding. In ^'^F 
our calculations yield an inversion of the Si/2 and ^5/2 
states, giving a 1/2 state as the ground state of ^^F. In 
order to gauge the effects of omitted short-range 3NFs, 
we vary the resolution scale A by a similarity renormal- 
ization group (SRG) transformation 19]. 

Figure [T7] shows the dependence of the neutron ^3/2 
single-particle state in ^^O as a function of the SRG cutoff 
A. In these calculations we used N + 1 = 15 oscillator 
shells with a fixed oscillator frequency oi fku — 26MeV. 
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FIG. 17: (Color online) Single-particle Energy (Eq. dSTj) of 
the neutron 0(3/2 state in ^*0 as a function of the resolution 
scale A of the SRG evolved chiral NN interaction. The calcu- 
lations employ A^ = 14 shells at fixed hu — 26 MeV. 



The (i3/2 single-particle state in ^''O depends weakly 
on the cutoff and stays bound for all choices of the cut- 
off A. This suggests that long-ranged 3NFs are needed 
to yield an unbound state and to place the drip line at 
^^O. Furthermore, the effect of the scattering continuum 
is expected to add some additional binding to this state 
(see Ref. [131 )• Thus, a realistic description of ^^O will 
result from a fine interplay between the scattering con- 
tinuum and 3NFs. Figure [T5] shows the cutoff depen- 
dence of the CCSD and A-CCSD(T) ground-state ener- 
gies for ^*0 (using the same model space and oscillator 
frequency as in Fig. [17)) . Note that there is no cutoff A 
that simultaneously reproduces the experimental binding 
energy of ^^O and the resonance energy of ^^O ground 
state. Figure fT9l shows that proton 51/2,^3/2, and d^/2 
single-particle states in ^"'O exhibit a considerable depen- 
dence on the cutoff A. No single cutoff reproduces simul- 
taneously the ground-state binding energies of ^''O, ^^O, 
and ^^F, again pointing to the importance of the omitted 
3NFs and to subtleties in its isospin dependence. 

In conclusion, chiral NN interactions alone do not re- 
produce the evolution of single-particle energies within 
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FIG. 18: (Color online) Ground-state energy of O within the 
CCSD (circles) and the A-CCSD(T) approximation (squares) 
as a function of the resolution scale A of the SRC evolved 
chiral NN interaction. The calculations employ in A^ + 1 = 15 
shells at fixed ^ = 26 MeV. 
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VI. RESOLUTION-SCALE DEPENDENCE FOR 
*^CA AND POWER COUNTING 



In this section we study the ground-state properties of 
^"^Ca using renormalized nucleon-nucleon interactions de- 
rived with the SRG method [T^ . The SRG method drives 
the Hamiltonian to a band diagonal form, and there- 
fore decouples low-momentum degrees of freedom from 
high-momentum degrees of freedom. The cutoff param- 
eter A determines the decoupling and sets the resolution 
scale or energy scale with which we can probe the struc- 
ture of a particular nucleus. However, this procedure 
induces three-body forces and forces of higher rank. At 
the two-body level, these interactions typically overbind 
medium-mass nuclei considerably, and show a strong de- 
pendence on the cutoff A. The dependence on the cutoff 
A on calculated observables gives an indication of the 
missing physics and on the rol e of many-body forces not 
included in the calculation; only the sum of all forces in- 
duced by the renormalization is truly independent of the 
scale of resolution or the cutoff. 

For the ground-state calculations of "'"Ca, we use SRG 
interactions at the resolution scales A = 2.5, 2.2, 1.9 fm~^, 
and compute the binding energy within the CCSD ap- 
proximation. We solve the CCSD equations in a spheri- 
cal harmonic oscillator basis. The SRG interactions are 
soft and the corrections due to triples clusters are found 
to be rather small. For ^"^Ca and a low- momentum in- 
teraction Viowfc with cutoff A = 1.9fm~^, for instance, 
the CCSD binding energy is E/A = 12.28 MeV and the 
CCSD(T) approximations yields an additional 0.29 MeV 
per nucleon |70j . Therefore, we limit the computations 
for the SRG interactions to the CCSD approximation. 
Our focus is on the saturation and convergence proper- 
ties of the SRG interactions in medium-mass nuclei and 
not on precision results. 



FIG. 19: (Color online) Single-particle energies (Eq. (|37|) ') of 
the proton 81/2,^3/2, and ^5/2 states in ^^O as a function of 
the resolution scale A of the SRG evolved chiral NN inter- 
action. The calculations employ A'^ -I- 1 = 15 shells at fixed 
&j = 26 MeV. 



the sd shell. Close to ^^O, the single-particle states are 
in semi-quantitative agreement with data, but missing 
3NFs become increasingly more important for neutron- 
rich isotopes of oxygen. It is expected, that a realistic 
description of single-particle states and the evolution of 
shell structure will require a theory which allows for a 
consistent and systematic inclusion of many-body corre- 
lations, 3NFs, and coupling to the scattering continuum. 
The beauty of coupled-cluster theory is that it allows for 
inclusion of all these ingredients in a simple and transpar- 
ent way, and we aim to investigate and predict properties 
of nuclei from the valley of stability to the very limits of 
matter taking all these ingredients into account in the 
near future. 



-610 




FIG. 20: (Color online) CCSD ground-state calculation of 
^"Ca for increasing number of oscillator shells, A'^ = 2n + Z, at 
fixed hw = 26MeV, using SRG evolved chiral NN interaction 
at resolution scale A = 1.9fm~ . 

Figures [201 Ell and [52] show the convergence of the 
ground-state energy of '^"Ca using SRG evolved interac- 
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ing the different SRG interactions. To examine the power 
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FIG. 21: (Color online) Same caption as in Fig. 1201 except 
that resolution scale is A = 2.2fm~^. 
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FIG. 22: (Color online) Same caption as in Fig. [20] except 
that resolution scale is A = 2.5fm~^. 



tions with resolution scales A = 1.9,2.2,2.5 fm^^, re- 
spectively. The ground-state energies computed for the 
various cutoffs are well converged in model space sizes 
comprising up to 13 major oscillator shells. As the cutoff 
A is lowered, these SRG interactions become increasingly 
soft, and we observe a faster convergence with increas- 
ing basis size. The computed ground-state energies show 
that the SRG NN interactions alone overbind consider- 
ably with respect to experiment. The results also exhibit 
a strong dependence on the resolution scale A, and the 
overbinding with respect to the experimental binding en- 
ergy of ""^Ca, which is approximately 342 MeV, increases 
with decreasing momentum scale. This strong resolution 
dependence is, of course, due to our omission of three- 
nucleon forces 3NFs (and more complicated many-body 
forces) and higher-order terms |66l. |67|. Likewise, with 
increasingly softer NN interactions, the proper reproduc- 
tion of saturation properties must come from the repul- 
sive character of the 3NF [64I . 

Table |ll] shows the binding energy per particle E/A, 
interaction energy per particle V/A, the momentum ex- 
pectation value Q, and the difference in binding energy 
per particle compared to experiment AE/A for '^"Ca us- 
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TABLE II: CCSD results for *"Ca with SRG interactions 
evolved from the chiral N LO nucleon-nucleon interaction to 
the momentum scale A. The binding energy per nucleon, and 
interaction energy per nucleon are denoted as E/A and V/A, 
respectively. Q denotes the expectation value of the momen- 
tum probed in this nucleus, while AE denotes the difference 
to the experimental binding energy. Energies are in units of 
MeV and momenta in units of fm~^. 

counting, we compute the expectation value of the inter- 
action energy via the Hellmann-Feynman theorem, and 
also deduce the average momentum Q from the expec- 
tation value of the kinetic energy. The results are given 
in Table |lTl The CCSD results are very well converged 
for momentum scales A = 1.9,2.2,2.5 fm^^. As already 
pointed out, our results are far from complete. Due to 
the restriction to NN forces, we are missing contributions 
of 3NFs at order N^LO, and NN, 3NF, and 4NF forces 
from higher order. The missing contributions from clus- 
ters due to the CCSD approximation are much smaller 
than the contributions related to neglected terms of the 
interaction. Chiral effective field theory puts the contri- 
butions of 3NFs at order {Q/\f. The SRG transforma- 
tion shifts high-momentum NN contributions to forces 
of higher rank but is not expected to destroy the power 
counting. Are our results consistent with this expecta- 
tion? To answer this question, we compute the ratio 
\AE/V\ / {\/Q)f for the ''°Ca nucleus and present the 
result in Table |TT1 The ratios are of natural size and sug- 
gest that the power counting is not violated in medium- 
mass nuclei. 

Let us discuss the implications of this finding. From 
a theoretical point of view, it is reassuring and satis- 
fying that the power-counting arguments seem applica- 
ble in medium-mass nuclei. From a practical point of 
view, this raises the question about the order necessary 
to achieve a desired precision. Low-momentum and SRG 
versions of NN interactions are quite popular since they 
facilitate considerably the solution of the nuclear many- 
body problem. However, our results also suggest that 
the momentum expansion in Q/A has to be driven con- 
siderably beyond the present order [Q/XY in order to ob- 
tain a small error estimate due to missing contributions 
from the nuclear interactions. Recall that the 4NF en- 
ters at order N'^LO. Excluding cancellation among 4NFs 
(which might perhaps be achieved through some fine- 
tuning of low-energy constants), the technical advantage 
of soft NN interactions might interaction wise be offset 
due to the technical difficulties that come with the inclu- 
sion of forces of increasingly higher rank. This potential 
problem has not been noticed in light nuclei. In light 
nuclei, the contributions of 3NFs are reduced at rather 
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small cutoffs due to fortuitous cancellations that are not 
expected from the power counting. Here, the lower cut- 
off benefits the handling of the NN force and it reduces 
the contributions of 3NFs. In light systems, there also 
seems to be only very little need for forces of rank higher 
than three. Table HIl hints that this might be different in 
'^^Ca. Here, and in agreement with the power counting, 
the contributions of 3NFs increase with the decreasing 
momentum scale of the underlying interaction. 

The results reported in this section suggest that there 
are at least two possible routes to ab-initio calculations in 
medium-mass nuclei. First, one can employ a low resolu- 
tion scale A which facilitates the solution of the nuclear 
many-body problem. The price tag to be paid might 
consist of a slower convergence with respect to the or- 
der kp/X, i.e., the evolution of four-body forces might 
be needed. Second, one can employ a high cut-off scale 
A^ and work on the technical difficulties related to the 
solution of the nuclear many-body problem. In this case, 
the full N'^LO interaction including 3NF forces might al- 
ready suffice to achieve accuracies of the order of 0.3 MeV 
per nucleon, and fine-tuning of poorly constrained coeffi- 
cients of the 3NF might further increase the accuracy. In 
following the second path, it might be attractive to em- 
ploy other techniques to tame the "hard" NN interaction. 
In the next section, we therefore explore the convergence 
properties of the G- matrix computed from the chiral NN 
interaction at order N'^LO. 



VII. CONVERGENCE PROPERTIES OF THE 
G-MATRIX 

In this section we calculate CCSD ground states of 
^He and ^^O using the in-medium G-matrix [22|. The 
G-matrix is an effective in-medium interaction that is 
computed starting from the bare nucleon-nucleon (NN) 
interaction using a Green's function approach with un- 
perturbed propagators [23]. It depends therefore on the 
starting energy ZU. The starting energy defines the energy 
of the incoming and outgoing single-particle states, and is 
normally set equal to the sum of the unperturbed energies 
of the interacting single-particle states. The G-matrix is 
defined in a model space P by summing ladder diagrams 
to infinite order, where the intermediate particle states 
are defined in the Q-space. Each ladder diagram scat- 
ters two particles from the P-space to the complement 
Q-space and back to the P-space. The number of in- 
teraction vertices in a ladder diagram gives the number 
of times the particles rescatter within the Q-space. It is 
clear that the G-matrix depends inherently on the model 
space and the interaction. The purpose of the G— matrix 
is to tame the hard core of the "bare" interaction mak- 
ing it suitable for many-body perturbation theory. In 
this section we investigate how the G-matrix behaves in 
ab-initio coupled-cluster ground-state calculations of "^He 
and ^^O as the model space P increases and the starting 
energy is varied over a wide range. 



In constructing the G-matrix, we start from the "bare" 
N^LO interaction (A^ == 500 MeVc"^) by Entem and 
Machleidt [Ij, \1^ and a spherical harmonic oscillator 
basis. In Figs. [23] and [Ml we show the CCSD results for 
^He and ^^O, respectively, starting from the calculated 
G-matrix. In the calculations we varied the starting en- 
ergy ID from -140 MeV to -5 MeV and studied the effect 
of this variation on the CCSD ground state energy. We 
increased the model space from 10 to 14 major oscilla- 
tor shells while the oscillator frequency was held fixed at 
hu = 20 MeV. The figures show that the dependence on 
the starting energy is reduced as we increase the model 
space. In the largest model space considered here, the 
starting-energy dependence is very mild. As ZU is var- 
ied over a range of 135 MeV, the binding energy of ^He 
changes only by about 0.25 MeV and of ^^O by about 
5 MeV. Furthermore, we see that at a fixed starting en- 
ergy, the binding energy converges rather slowly with in- 
creasing size of the model space. As the model space 
P increases the complement space Q in which the par- 
ticles scatter is reduced, and the G-matrix should there- 
fore converge to the underlying "bare" nucleon-nucleon 
interaction in an infinite model space. This expectation 
is compatible with our CCSD results for the ground- 
state energies of ''He (about —24 MeV) and ^^O (about 
— 108 MeV) using the "bare" chiral NN interaction (see, 
e.g.. Figs. [T] and ^. Figures [^ and [Ml show that the 
ground-state energy of the G-matrix converges from be- 
low as the model space is inceased in size, and that the re- 
sult of the "bare" interaction is approached in very large 
model spaces. 
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FIG. 23: (Color online) CCSD binding energy for *He using 
a G-matrix and starting energies Zj in the range from —140 
MeV to —5 MeV. The oscillator spacing is held fixed at hu — 
20 MeV. 

In conclusion, we observe slow convergence for the 
CCSD ground-state energy with increasing model space 
with the G-matrix defined from the "bare" N'^LO interac- 
tion. The dependence on the starting energy disappears 
in sufficiently large model spaces, and the coupled-cluster 
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FIG. 24: (Color online) CCSD binding energy for ^"O using 
a G- matrix and starting energies U in the range from —140 
MeV to —5 MeV. The oscillator spacing is held fixed at fia; = 
20 MeV. 



results with the G-matrix converge towards coupled- 
cluster results for the free nucleon-nucleon interaction. 
In the case of the N'^LO interaction, it is clearly bet- 
ter to start from the free interaction. It remains to 
be seen in the case of harder interactions, like the Ar- 
gonne V18 interaction, whether convergence is better us- 
ing the in-medium G-matrix as a starting point for ab- 
initio coupled-cluster calculations. 

Let us briefly contrast the G matrix to low-momentum 
NN interactions a la Viow k or from the similarity renor- 
malization group. The latter two preserve phase shifts 
from NN scattering and bound-state properties up to the 
chosen momentum cutoff. They soften the NN interac- 
tion and - when not augmented by the induced 3NF - 
tend to considerably overbind heavier nuclei at lower cut- 
offs around A w 2 fm~^ or so. In a finite model space, the 
G-matrix always exhibits a dependence on the starting 
energy, it converges slowly to the "bare" interaction, and 
it somewhat overbinds the nuclear many-body system 
compared to the "bare" interactions that is used in its 
construction. The authors speculate that this last point 
might seem attractive with a view on practical applica- 
tions, as "bare" nucleon-nucleon interaction models such 
as the Argonne interaction, the CD-Bonn interaction and 
the chiral interaction [l3| somewhat underbind nuclei. In 
other words, the slight overbinding of the G-matrix often 
results into very reasonable energies for heavier nuclei, 
and the G-matrix simply works very well in practical ap- 
plications. 



VIII. CONVERGENCE PROPERTIES OF Vucom 

Let us also consider the saturation and conver- 
gence properties of the Vucom nucleon-nucleon interac- 
tion j24| . The interaction Vucom is obtained from a uni- 



tary correlation operator method (UCOM) that softens 
the short-range repulsion of the initial "bare" interaction 
by a similarity transformation with a unitary operator, 
which is explicitly designed to remove the hard-core and 
the short-ranged tensor components of phenomenological 
NN interaction. In this section we consider the UCOM 
interaction obtained from the "bare" Argonne VIS inter- 
action. 

The Vucom interaction was initially applied in mean- 
field methods. These applications were very successful. 
Very reasonable binding energies and saturation proper- 
ties could already be obtained within the Hartree-Fock 
approximation. Only later did this interaction see ap- 
plications that included the inclusion of many-body cor- 
relations within the RPA, many-body perturbation the- 
ory [7l| , the coupled-cluster method |72|] , and within the 
no-core shell model |73l|. (For a recent review, we refer 
the reader to Ref. [TJI-) Those later applications suggest 
that the Vucom interaction is soft yet with appealing 
saturation properties, and without the need for sizeable 
3NFs. This situation seems surprising. On the one hand, 
the Vucom interaction is similar in its technical construc- 
tion to the SRG interactions jT^- On the other hand, 
the Vucom interaction seems to differ significantly in its 
saturation property from the SRG nucleon-nucleon inter- 
actions which tend to overbind nuclei at low momentum 
cutoffs. We address this puzzle in what follows, and per- 
form structure calculations with the UCOM interaction 
in large model spaces and for various oscillator frequen- 
cies. 

We employ the Vucom interaction and compute the 
ground states of ^^O and ^"Ca within the CCSD ap- 
proximation. Our calculations employ very large model 
spaces and a considerable range of oscillator freq uencies, 
and thereby differ from the previous studies |7ll - l73{ : see 
also Ref. f25| . Figures [25] and [26] show the ground-state 
energy of ^^O and ^°Ca, respectively, as a function of 
the employed oscillator frequency uj for different sizes of 
the model space {N + 1 denotes the number of oscillator 
shells.). The convergence with respect to increasing size 
of the model space is slow, and we are unable to achieve 
convergence in up to 14 major oscillator shells. The 
results obtained in the largest model spaces also show 
that the Vucom interaction considerably overbinds ""^Ca. 
Within the CCSD approximation, we find a ground-state 
energy of about —133 MeV for -'^^O and of about —400 
MeV for ^°Ca. The inclusion of triples corrections within 
the A-CCSD(T) approximation yields a ground-state en- 
ergy of -140.99 MeV for ^'^O in 13 oscillator shells and 
an oscillator frequency hoj — 42 MeV. 

We checked our results as follows. First, our CCSD 
results agree with those reported for '*He in Ref. [72| 
within about 50 keV in model spaces of up to eight oscil- 
lator shells. Second, we also computed the ground-state 
energy of *°Ca in third order many-body perturbation 
theory. The relevant diagrams which are included can 
be found in Ref. [33] . Our results are shown in Fig. [57] 
as function of the oscillator energy hu and the number 



19 



M 



-105 



-110 



-115 



-120 



■125 



-130 





1 1 1 1 1 1 


1 1 1 


1 1 


; 






•-•N = 7 




: 


- 




■■■■N = 8 




- 




16 


♦■-♦N-9 




■ 




withV^,^^ 


^iN= 10 




: 


_ 




<-<N= 11 




_ 




•\^ 


N=12 




: 




^\^ 


► >N=13 




: 


T 




— Exp. 


__ 


- 




- 


^-:-^:■■■■■ 




...■■■■ 


- 


■ ■■■ 






~--. ^~ --♦-- 






■ 




■'•■-- ~~-*-_ - 


♦ »_ 


^ 


: 


- 


''^■M-. ~" ~ 


i A- 


_ -A 


- 




---:fc:-:;:t 




: 


> 







1 1 1 


1 1 


■ 



"^■^^0 22 24 26 28 30 32 34 36 38 40 42 44 
fico (MeV) 

FIG. 25: (Color online) CCSD ground state energy of ^"0 
with VucoM as a function of the oscillator spacing huj and 
the size of the model space A''. The experimental ground- 
state energy is given by the solid line. 
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FIG. 27: (Color online) Ground-state energy for ^"Ca using 
many-body perturbation theory to third order in the interac- 
tion VucoM- The results are presented as a function of the 
oscillator spacing &j and the size of the model space A''. As in 
the CCSD calculations, a Hartree-Fock basis has been used. 



This speculation is borne out by two observations. First, 
the comparison of the UCOM interaction with the SRG 
interaction in Fig. 8 of Ref. [75^ shows that the UCOM 
interaction does not decouple well in an oscillator basis. 
Second, the energy minima in Figs . B5l [^ and \T7\ shift to 
higher frequencies hut as the size of the model space is in- 
creased, as one would expect based on the estimate (|36)) 
for an interaction with a large momentum cutoff. Our re- 
sults also suggest that - similar to other SRG interactions 
and low-momentum interactions - the UCOM interaction 
has to be augmented by sizeable many-body forces [7q . 
The need for these forces is not seen if one restricts the 
calculations to fixed model spaces where model space pa- 
rameters such as the oscillator frequency and the number 
of shells introduce additional momentum cutoffs. 



IX. SUMMARY 



of shells N. The results from many body perturbation 
theory confirm the trend seen in the coupled-cluster cal- 
culations, i.e., they exhibit a rather slow convergence in 
terms of the number of shells N and overbind the nu- 
cleus '"'Ca. Reference (tH also solved the UCOM po- 
tential within many-body perturbation theory, albeit for 
oscillator frequencies below /zw « 20 MeV or so. In this 
regime, the binding energy is closer to the experimental 
value as the small oscillator frequency acts as a cutoff in 
momentum space. 

We speculate that the Vucom interaction still exhibits 
relatively long tails in momentum space, which prevent 
a decoupling of low- and high-momentum modes, or that 
peculiarities of the corresponding cutoff yield a slow con- 
vergence of the UCOM interaction in an oscillator basis. 



We presented a spherical formulation of the coupled- 
cluster method for the computation of energy spectra in 
nuclei with closed subshells and their neighbors. This 
method has been used to solve nuclear structure prob- 
lems in model spaces consisting of up to 20 oscillator 
shells, which allows us to obtain well-converged results 
for nucleon-nucleon interactions derived from chiral ef- 
fective field theory. We find that chiral NN interactions 
saturate nuclei such as ^^O and ^^Ca within 0.5 MeV per 
nucleon compared to data. The two-particle-two-hole 
clusters provide about 90% of correlation energy, with 
the approximation of triples clusters accounting for the 
remaining 10%. We investigated the shell evolution in 
neutron-rich isotopes of oxygen and fluorine and found 
that nucleon-nucleon interactions alone fail to describe 
the experimentally observed (sub) shell structure. 
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We also employed similarity renormalization group 
transformations of a "bare" chiral interaction. At the 
considered resolution scale, we obtained overbinding of 
up to several MeV per nucleon. Again, the missing con- 
tributions (when compared to experiment) are of natural 
size in the power counting of effective field theory. These 
results suggest that the systematic approach to nuclear 
structure can (in principle) be extended to medium-mass 
nuclei. We studied the G-matrix approach and found 
a very weak starting-energy dependence in large model 
spaces combined with a rather slow convergence with re- 
spect to increases in the size of the model space. For 
the UCOM interaction, we find a slow convergence with 
respect to the size of the model space and an overbind- 
ing for *'^Ca that is comparable to other low-momentum 
interactions. 

We presented an in-depth study of the center-of-mass 
problem, and demonstrated that the wave functions of 
the intrinsic Hamiltonian factorize to a very good approx- 
imation into an intrinsic wave function and a Gaussian for 
the center of mass. While an analytical understanding of 
this behavior is still lacking, our calculations for ground 
and excited states and for a virtually exactly solvable 
toy model indicate that calculations in sufiiciently large 
model spaces - when based on the intrinsic Hamiltonian 
- do not suffer from a center-of-mass problem. 



Appendix 



The two-body amplitudes are diagonal in the total an- 
gular momentum, that is, we have Jij = Jab- The fi- 
nal matrix elements are independent of M. The labels 
a,b,c,d . . . refer to particle states while i,j,k,l,... are 
hole states. 




FIG. 28: (Color online) Diagram contributing to the t1j{pp) 
amplitudes in particle-particle coupled scheme. 



In Fig. [28] we give a diagram contributing to the T2 
equation, and the corresponding algebraic expression in 
the J— coupled scheme is (we list only the angular mo- 
mentum related part), 



X. COUPLED-CLUSTER DIAGRAMS IN A 
J-COUPLED SCHEME 



lab 



(Jij) ^^{{Jajb)Jab\V\{jcjd)Jab){{jcjd)Jab\t\{jtjj)Jab) 



cd 



In this section we present expressions for several dia- 
grams involved in our computation of the CCSD or A- 
CCSD(T) equations that involve various one-body, two- 
body, and three-body amplitudes and operators. 



Examples of CCSD diagrams in a J-coupled 
scheme 



We employ the following shorthands for the Ti and 
T2 amplitudes, namely tf and tfj, respectively. To the 
latter we add the total two-particle angular momentum 
Jij in order to indicate a two-body wave operator with 
hole states ij coupled to a two-body angular momentum 
Jij, that is 



t^jU^j) ^ ((j-3b)J^MJ^Jl)J^: 



E 

mamhrriiraj 






The above expression is obtained by identifying the dia- 
gram as a ladder cutting apart the lines connecting the 
various vertices. This technique is described in detail in 
Ref . [73] . Several applications of this technique are listed 
also in Ref. [22]. This is the most expensive diagram in 
the CCSD approximation and scales as n^n^, where Uo 
is the number of occupied j'— orbitals and n„ is the num- 
ber of unoccupied j-orbitals. As seen, the diagram can 
be calculated by a matrix-matrix multiplication utilizing 
efficient BLAS routines [Zlj. The diagram in Fig. B^ also 




{{iama){jbmb)\t\{jim.i){ijmj)) , 



(38) 



FIG. 29: (Color online) Diagram contributing to the tij (ph) 
amplitudes in particle-hole coupled scheme. 



with mi,mj, etc. being the magnetic quantum num- 
bers of the corresponding single-particle angular mo- 
menta ji,jj, etc. The coefficients C are the standard 
Clebsch-Gordan coefficients. 



contributes to the T2 equation, and is easiest calculated 
in a particle-hole coupled scheme so that intermediate 
sums of 9j— symbols are avoided. The transformation of 
the tfHJij) amplitudes to a particle-hole coupled scheme 
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is given by 

{UaJ~')Ja^\t\{JJJ^')Ja^) ^ 

{{jajb)Jab\t\{jijj)Jab) , (39) 

and the transformation back to particle-particle coupling 



[4] 



{{jajb)Jab\t\{jijj)Jab) = 

^(2J„ + l)(_i)^-»+i.+.'-+.^=. 



Ji Jj 'Jab 
Jb Ja 'J ai 



(40) 



The diagram in Fig. [29] can now easily be calculated in a 
particle-hole coupled scheme as, 

{{Jaj-^)Ja^\i\{3,h^)Ja^) ^ 



ck 



{{]kJ-^)Ja^\V\{J,J^^)Ja 



(41) 



where again an efBcient BLAS [78| matrix-matrix multi- 
plication routine can be used. 

In Fig. 1301 we give a diagram contributing to the Ti 
equation. The algebraic expression of Fig. [30] in a J- 




FIG. 30: (Color online) Diagram contributing to the t^ am- 
plitudes. 

coupled scheme is, 

*"^ESTE^"(J..)x^ 

Jik cfc 

where x^ is a one-body operator. 



B. Examples of A-CCSD(T) diagrams in a 
J-coupled scheme 

Figure [ST] shows an example of diagrams with a three- 
body wave operator that give rise to contributions to 
the final energy. These diagrams arise from the so-called 
A-CCSD(T) approximation discussed in this work. The 
lower part of these diagrams is given by the three-body t- 
amplitude labeled t^^^. It is antisymmetrized and in the 




E 



FIG. 31: (Color online) Diagrams for the A-CCSD(T) approx- 
imation. 



A-CCSD(T) approximation it is represented by a con- 
nected three-body contribution consisting of a two-body 
amplitude tfj and a two-body interaction vertex plus a 
two-body amplitude part if.* times a one-body amplitude 
f^ . These two contributions are depicted by the diagrams 
to the right of the equality sign in Fig. [32] The three- 
body amplitude t^^^ is shown to the left of the equality 
sign in Fig. [32] 



.V 



FIG. 32: (Color online) Diagrams for the A-CCSD(T) approx- 
imation. 

A diagram like the first to the right of the equality sign 
in Fig. [35] can be redrawn as shown in Fig. 1331 Here we 



FIG. 33: (Color online) Connected part of three-particle- 
three-hole diagrams. 

have stretched the hole lines in the diagram to the right 
of the equality sign. Such a diagram with, for example, 
an intermediate particle state can easily be calculated in 
an angular momentum coupled scheme. The expressions 
can be derived using the methods discussed in Ref. (77| . 
All three-body diagrams like those shown in Fig. [33] have 
been listed in Ref. [79} . 

In order to derive these expressions, one needs to spec- 
ify the given coupling order for the angular momenta. 
Here we choose to couple our angular momenta as 

|([ja -^ 3b]Jab --^ ic)J) = 

E Uarriajbrnb \ JabMab){JabMab3cTnc \ J M) 

marrif^mc 

iJarUa) «) iJbTOfc) (K) iJcTOc) , (42) 

and indicated in Fig. [34l 

Note that the two-body intermediate state is antisym- 
metric but not normalized, that is, the state which in- 
volves the quantum numbers ja and jb- We will hereafter 
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3b 



3c 



Jab 



3d 









FIG. 34: The coupling order ([> -> jb]Jab -> jc)J- 



assume that the intermediate two-body state is antisym- 
metric. With this couphng order, we can rewrite the 
general three-particle Slater determinant as 



<^{1,2,3) ^ A\{[ja ^ 3b]Jab ^ 3c)J), 



(43) 



with an implicit smn over Jab- The final Slater determi- 
nant is 



Ji 



3j 



3k 



FIG. 35: Three-body diagram with two-body vertex and am- 
plitude and a particle state as intermediate single-particle 
states. The sum over d runs over all possible quantum num- 
bers of the intermediate particle(hole) state. The algebraic 
expression of the above diagram, in an angular momentum 
coupled scheme, is given in Eq. (1441) . 



^I([ja -> jb]Jab -^ 3c)J) = 
^l([ia -^ 3b] Jab -^ jc)J) 



J2i-iy^+'^ + '''^ + ''-JabJac 



J2{-iy'^JabJbc 



./be 



ja jb Jab 
'J Jc ^bc 



I ja jb Jab I . 
I '-' Jc ^ac J 

lib a -> jb]Jab -> jc)J) 



Again, we employ the shorthand tf^(Jij) to indicate a 
two-body wave operator with hole states ij coupled to 
a two-body angular momentum Jij, that is, i^^j(^ij) = 
{Ua3b)Jij\t\{jijj)Jij)- The diagram with a particle inter- 
mediate state is shown in Fig. |3S] with its angular mo- 
mentum representation, 



Here, we used the shorthand J = V2J+T. With this 
coupling order, one can then compute all three-body di- 
agrams like those listed in Fig. |33l 

Here we give an example of the closed form expres- 
sion for the angular momentum recoupled part of a se- 
lected three-body diagram with two-body amplitudes tfj. 
These are diagrams which start with a three-body state 
but have two-body wave operators and end in a con- 
tribution to a three-body state with a two-body ver- 
tex after the two-body amplitude, as shown in Fig. 1351 



^^(_ip+.c+.. + J..J^^J^^.j2^, 



Jbc d 



Ja Jb 'Jab 
Jc J 'Jbc 



Ja Jd 'Jij 
jk J Jbc 



{{jbJc)Jbc\V\{MJk)JbcKMj)-i^^) 



The expression for this diagram is obtained by opening 
up the intermediate three-particle state jajdjk and re- 
coupling the angular momenta jb and jc to yield a final 
two-particle angular momentum J^c- This applies also to 
the single-particle angular momenta jd and jk ■ These two 
single-particle angular momenta couple to the same final 
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two-particle angular momentuni Ji^c- These recouplings 
are reflected in the two 6j-synibols in the above expres- 
sion. The expression for the corresponding diagram of 
Fig. [35] with a hole intermediate state is obtained by re- 
placing the particle labeling jd with j;. The subscript 
/ refers to the fact that this is a single-hole state. The 
internal hole line, however, gives rise to a factor —1. 

In total there are, due to the antisymmetry of the 
three-body wave function, nine diagrams with an inter- 
mediate particle state and nine diagrams with an inter- 
mediate single-hole state. The angular momentum ex- 
pressions for these diagrams are listed in Ref. [73] ■ The 



expressions for diagrams like the rightmost one in Fig. [32] 
yield similar expressions and are easy to compute. They 
represent a two-body part coupled with a one-body part. 
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